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1 Introduction 



Consider a hexagonal lattice with M vertices. In how many ways is it possible to 
colour the links of the lattice with three colours, A, B and C so that no two links 
which meet at the same vertex carry the same colour. This is the classical three-colour 
problem solved by Baxter in 1969 [|I|. The problem is equivalent to a loop gas problem, 
namely the problem of enumerating the number of ways one can cover the lattice with 
closed, self-avoiding and non-intersecting loops which come in two different flavours, 
demanding that each vertex of the lattice belongs to a loop. This follows by noting 
that sequences of, say, B- and C-coloured links form loops which exactly have the 
properties described above The loop gas model, on the other hand, is a special 
version of the 0(2) model 0. Using the equivalence with the loop gas model Baxter 
furthermore proved that his problem was equivalent to the problem of colouring the 
faces of the hexagonal lattice with four different colours so that adjacent faces have 
different colours The dual lattice of a hexagonal lattice is a triangular lattice. 
Baxter's three-colour problem can hence also be formulated as the problem of counting 
the number of ways of colouring the links of the triangular lattice so that the three 
sides of any triangle have three different colours. In addition, it has recently been 
shown that Baxter's three-colour problem is equivalent to the problem of counting the 
different foldings of the regular triangular lattice |^. 

In this letter we study the generalisation of Baxter's three-colour problem to a 
random triangular lattice. This can be viewed as coupling a hitherto unexplored type 
of matter to two-dimensional quantum gravity. Whereas the regular triangular lattice 
has six triangles meeting at each vertex, on a random triangular lattice any number 
of triangles can meet at a given vertex. However, the random lattice three-colour 
problem shares many features with its regular lattice version. For instance,the random 
lattice three-colour problem is equivalent to a loop gas model, namely a version of 
the 0(2) model on a random lattice where only loops of even length are allowed. 
Furthermore, as in the regular lattice case, the random lattice three-colour problem is 
equivalent to the problem of counting the possible four- colourings of the faces of the 
dual lattice. However, the random lattice three-colour problem can not be given an 
interpretation as a folding problem. This requires that only configurations where an 
even number of triangles meet at a given vertex are allowed. Finally, let us mention 
that the three-colour problem on a random lattice has recently attracted attention as 
a means of describing vertex models on random graphs 
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2 The Model 



The possible three-colourings of the random triangular lattice are generated by the 
following matrix model 



Z{g) = / dAdBdC exp 

Jnxn 



-NTt ( - 
,2 



A'' + B^ + C^ [ABC + BAC] (2.1) 



where A, B and C are hermitian NxN matrices. More precisely we have the following 
expansion of the free energy, F{g) = log Z{g) 



F{g) = Y: N-''F,{g), F,{g) = ^ A4(2t;) g^ 



(2.2) 



h=0 



where N'h{2v) is the number of closed, connected, three-coloured random triangulations 
of genus h consisting of 2v triangles. (Here it is understood that a given triangulation, 
T, is counted with the weight l/|Aut(T)| where |Aut(T)| is the order of its automor- 
phism group.) It is important to note that it is necessary to take into account both of 
the terms Tr (ABC) and Tr (ACB) in ( |2.1| ) if one wishes to allow for any cyclic order 
of the colours around a given triangle. If one considers only triangulations where the 
colours occur in the same cyclic order for all triangles only one term is needed and the 
model simplifies drastically. The model with only one term, Tr (ABC), was studied by 
Cicuta et. al. who extracted numerically the value of the critical coupling constant and 
the critical index, ■jstr 0- Later, it was shown that the model could be mapped onto 
an 0(1) model on a random lattice with a non-polynomial potential and this led to an 
exact solution of the model showing that the critical indices of the model coincide 
with those of the Ising model on a random lattice. 

In stead of the model (p.l|) we shall consider the more general model 



NxN 



dA n dB^ dCi 



i=l 



{1 1 " " 11 

-iVTr -V{A) + -Y,[Bf + C^)-Y.{ABia + AQBi) (2.3) 
.9 i=l i=l J J 



This model is not more difficult to treat than the model ( |2.1| ) and reduces to the latter 
for n = 1, and V{A) = \A^ , i.e.Q 



Z{g) = i^f Z{g) 



n=l,V{A)=^A'2 ■ 



(2.4) 



^ We note that in analogy with the n = 1 case discussed above, for general n the model ( |2.3| ) with 
only one term X]r=i {ABid) term can be mapped onto an 0{n) model with a non-polynomial 
potential and solved exactly. In particular, the critical indices of the model can be shown to coincide 
with those of the 0{n) model [Q. 
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Integrating over the C-matrices in (| 2.3| ) we get 
Z{g)= [ dAUdBiexpl-NTi 

JNxN Z_i 



9 j=i i=i 



(2.5) 

In this letter we shall concentrate on solving the counting problem for triangulations 
of spherical topology. We will derive an equation which allows us to calculate, in 
an efficient way, the genus zero contribution to all correlators of the type (;^Tr A"') 
iteratively in g. It will be clear, however, how to extend the idea to obtain the higher 
genera contributions. We note that just knowing (;^Tr A^) is enough to solve the 
three-colour problem, since we have for n = 1, V{A) = ^A^ 

(4 Tr A') = %4z ^ogZ{g) = 2/-^ (i log(^) + log Z{g)\ (2.6) 



'N ' N'^ dg ^ " dg\2 ^"-^^ m 

and in particular for genus zero (cf. equation (|2.2|) ) 

1 oo 

(- Tr A')^=o = g + 2gJ2 ^o(2i;) v g^ = gT,. (2.7) 

v=0 

Baxter solved the regular lattice three-colour problem by a transfer matrix method. 
Recently a transfer matrix formalism for random triangulations has been invented |^ 
and it is natural to ask whether this formalism can be applied to the present model. It 
indeed can. A version of the transfer matrix formalism on random lattices, applicable 
to loop gas models, was presented in and can be applied to the three-colour problem 
exploiting its equivalence with a variant of the 0(2) model on a random lattice. In this 
approach no reference to any matrix model description is needed. However, we shall 
stay within the matrix model approach because this gives a faster way of deriving the 
equations we need. 

In the subsequent section we shall write down the saddle point equation correspond- 
ing to the matrix integral ( |2.5|) . This exposes the analyticity structure of the model 
and shows why the full three-colour problem is so much more complicated than the 
restricted one. Actually one can derive the equation we are ultimately after entirely by 
analyticity arguments, based on the saddle point equation but we shall evoke another 
line of action, namely the loop equation method. This method has the advantage that 
it is immediate to see how to proceed to higher genera. 

3 The saddle point equation 

Integrating over the B- matrices in (|2.5| ) we get 

Z{g)= [ dAe"^^^^^^ (det(l®l + l(giA + ^l)(lOl-l®A-^l))"^. (3.1) 

J NxN 
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Furthermore, integrating over the angular degrees of freedom leaves us with the fol- 
lowing integral over the eigenvalues of the matrix A 

N 

Z{9)= j X{^KX{{\,-\kfX{{{l + k + \rn){l-k-\rn)r'''e-TT..yi^^), (3.2) 
i=l j<k l,m 

In the limit N ^ oo the eigenvalue configuration is determined by the saddle point of 
the integral above. The corresponding saddle point equation reads 



2Y.^-nY. ^ ^ 

j^i -\ J L 1 + Ai + Xj 1 - Ai - Xj 



ly'(Ai) Vi = l...iV. (3.3) 



Following reference |]TD| we can introduce an eigenvalue density p{z) = ^{^ ~ ^i) 
which in the limit N ^ oo becomes a continuous function, and a resolvant W{z) = 
Si tta" ^ S dX As is clear from the expression (|3.2|) , the model becomes singular 
if one of the eigenvalues approaches —\ or \. The solution we are interested in corre- 
sponds to the situation where the eigenvalues live on only one interval [a, h] C [— |, |] 
or equivalently the situation where W{z) is analytic in the complex plane except for a 
cut [a, h] C [— I, |]- Written in terms of the resolvant the saddle point equation reads 

W{z + i{]) + W{z-iQ) + n{W{l-z) + W{-l-z)} =^V'{z)- ze[a,h]. (3.4) 

This equation describes how W{z) transforms when z crosses the cut and enters into 
another sheet. In the second sheet, W{z) is a combination oiW{z), W{1 — z) and 
W{—l — z). Thus in the second sheet there are three cuts: [a, 6], [—1 — 6, —1 — a] and 
[1 — 6, 1 — a] . Crossing again these cuts, we generate an increasing number of cuts in 
the next sheets, which means that W{z) is defined on a Riemann surface of infinite 
genus and with an infinite number of cuts in each sheet. As opposed to this, for the 
1-matrix model (n = 0) W{z) has only one cut and two sheets and for the 0{n) model 
on a random lattice W{z) has only two cuts in each sheet 0]. Here we need to consider 
and infinite series of cuts {Ik} given by 



a 



^^\h^% a^^^ = k + {~lfa, h'^^^ = k+{-lfh (3.5) 



and we note that the critical situation referred to above, where one of the eigenvalues 
approaches — ^ or i, corresponds to the situation where all the cuts merge. One might 
expect, in analogy with what was the case for the 0{n) model on a random lattice, that 
this type of critical behaviour can only be realized for a certain range of n- values 
We note that when this type of critical behaviour is realized the scaling behaviour of 
the eigenvalue distribution in the vicinity of the endpoints of its support will be as for 
the 0{n) model. For instance, if we consider p{z) or equivalently W{z) in the vicinity 
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of 2 = I we can view the term 1 — z) as being regular and the saddle point 
equation reduces to that of the 0{n) model on a random lattice. 

Let us now show how it is possible, using the saddle point equation, to build from 
W{z) a function which has no cut at all. First we define 

Zk = {-l)Hz-k) (3.6) 

Then we can write the saddle point equation as 

W{zk + lO) + W{zu -iQ)+n {W{zu+i) + W{zk^^)} = -V'{zk), z G K, h]. (3.7) 

9 

Multiplying this equation by W{zk + zO) — W{zk — iO) we find that the function 

g{z) = W\zk) + nW{zk) {W{zk+i) + W{zk-i)} - -W{zk)V\zk) (3.8) 

has no cut along (while it of course has cuts along l^-x and /fc+i). Now it is easy to 
see that the following function, S'(2;),has no cut at all 

OO 1 

S{z)= J2 W^iz,) + nWiz,)W{zk+,)--iWizk)V'izk)-Rizk)). (3.9) 

k=—oo 3 

Here R{z) is the polynomial part of W{z)V'{z) which we have subtracted in order to 
ensure that the sum converges. Since S{z) can have no singularities apart from the 
above mentioned cuts it must be analytic in the whole complex plane. Furthermore, it 
is easy to see that S{z) fulfils the following relations 

Sil-z)=S{z), S{z + 2) = S{z). (3.10) 

From the periodicity relation it follows that S{z) can be written as a Fourier series 

/•3/2 

Siz) =Y,Spe'''P', ^P = J/ S{z)-'''P'dz. (3.11) 

Now it is actually possible by pure analyticity arguments to show that Sp = Vp, i.e. 
that 5* (2;) is identically equal to zero. However, we shall derive this result by another 
method, namely by means of the loop equations of the model. These equation have 
the advantage that they contain also information about higher genera contributions. 

4 The loop equations 

The loop equations simply express the invariance of the matrix integral ( p.5| ) under 
analytic redefinitions of the integration variables. Let us introduce the notation B = 



5 



{Bi, . . . , Bn) and let us define 

where the subscript conn refers to the connected part. We note that the leading 
contribution to all of the above listed correlators is of the order N^. Now, let us 
perform the following redefinition of the field A 

A^A + e^. (4.4) 

This gives rise to the identity 

- {V'{z)W{z) - R{z)) = 2H{z) + 2zW2{z) - 2Tb + W^z) + -^W{z, z). (4.5) 

Furthermore, considering the following redefinition of the field B 

B ^ B + e^—-B ^— (4.6) 

z - A X - A 

we find 

H{z) + H{x) + (x + 2z)W2{z) + {z + 2x)W2{x) - 2Tb 

+ {l-{z + xf)F{z, x) - nW{z)W{x) - -^W(z, x) = 0. (4.7) 

We see that as announced the loop equations (as usual) permit a study of higher genera 
contributions. However, in the following discussion we shall restrict ourselves to genus 
zero. First, we note that equation ([4.7| ) simplifies considerably when z + x = ±1 = 6. 
Then, setting x = 6 — z and combining (|4.5| ) and (^]3) we find 

26 (W2(z) + W2(S - z)) = --{V'{z)W{z)-R{z) + V'{6~z)W{6-z)-R{6-z)) 

9 

+W\z) + 2nW{z)W{6 - z) + W'^{5 - z) (4.8) 

Replacing z hy Zk, and 5 by (—1)*^ we have 5 — Zk = Zk+i and summing over k (from 
ki to k2 — 1) we get 



k2-l 

E 

k=ki 



W\zk) + nW{zk)W{zk+i) - - {V'{zk)W{zk) - R{zu)) 

9 

\ (w\zk;) - ^ {y\z,,)w{zk,) - R{z,,))^ 

-\ (w^'(^fe) - ^ {V\z^,)W{z^,) - R{z,,)) 



+ {{-lf^W2{z,J - (-l)^W2(z,j) . (4.9) 
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The right hand side vanishes when ki — > — cxo and ^2 — > +00, while the left hand side 
is precisely the function S{z) of equation (|3.9|). Therefore we have 

S{z) = 0. (4.10) 

We thus have a closed (non-local) equation for W, which could in principle allow us to 
retrieve W{z). 



5 An equation for the moments 

We shall now write equation ( [4 .101 ) in a more manageable form, namely as an equation 
for the moments tm = (j^TrA^Y related to W{z) by 



m=0 ^ 



We shall use the normalisation condition to = 1 which corresponds to requiring that 
the eigenvalue distribution is normalised to one. 
First, we note that we can write 

00 

Siz) = Y.fizk) (5.1) 

fc=0 

with 

71/ 1 

f(z) = W\z) + -W(z)(W(l -z) + W(-l - z)) - -(V'(z)W(z) - Riz)). (5.2) 

2 g 

Next, we introduce the following functions 

U^)=E4 (5.3) 

k 

C (1 - ZkF T ^k (-1 - Zkf 

The functions Cm (-2) can be determined explicitly (by Fourier transform). One has 

Ci(^) = C2{z) = -4^, (5.5) 

smvrz sm irz 

1 d^ 

Cm+2{z) = —, r-^T^Crn{z). (5.6) 

m(m + 1) dz^ 

In other words 

/ TT \™ 

C,m{z) = [- P^(sin7r2;) (5.7) 

V sm TTZ / 
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where Pm is a polynomial of degree less than or equal to m — 1. Actually, we can give 
an explicit expression for this polynomial, namely Pm{s) consists of the first m — 1 
powers of s in the power series expansion of I — -. — for s in the vicinity of zero. The 

tr f y arcsm s J •' 

functions Cm,p(-2) and Cm,p(^) can be expressed in terms of the Cfc(-2) as follows 

WW = |:("\":i"')c.(.)+E("^+^:;-')a(.) (5.8) 

+ E ( : 1 " ^ ) {-^r^'^'UA (5.9) 

Now, starting from the expression (|5.1| ) and choosing the potential of our model as 

3 

we can write 



V{z) = Y.^z^ (5.10) 



n — if °° ^ 

k,l k,l 9 \ k=0 l=k J 



0. 



Using relations ( |5.8| ) and (p.9|), and identifying the coefficients of each C,rn we get for 
m > 1 

]^ degV m-2 oo 2fc+l / 2k + 1 \ 

- E 9ktm+k-l = E ^ktm-k-2 + ^'T- E E Oi- _L 1 7 ^m+«-1^2fc+l-i • (^-H) 
9 k=0 k=0 k=0 1=0 \ L-L J 

This equation is equivalent to equation ( [1.10| ) and contains all information about the 
full non-perturbative solution of the model. Lacking a means of solving this equation 
exactly we shall in the next section describe how a perturbative solution can be found. 

6 Perturbative solution 

Let us specialise to the case of a Gaussian potential, i.e. V'{z) = z, but keep n arbitrary. 
Setting n = 1 at any stage of the calculation then brings us to the three-colour problem. 
For a quadratic potential, obviously t2m+i = 0, Vm, and the function W{z) is an odd 
function of z which is analytic in the complex plane except for a cut of the type 
[—a, a] C [— |, |]. The relation ( p.ll[ ) in this case reduces to 



1 . ^ ^ I 2k + 21 + I 



t2m = y2 hkhim-k-l) + 2n E ( 07 I l' 1 ^2(m+fc)^2Z, 171 > 1. (6.1) 

9 k=o k,i=o \ ^'^^^ I 
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Setting t2m = g"^Tm we have 



m— 1 



Tm=Yl TkTm-k-1 + 2ngJ2{ J. ^J' g'+'T^+kTi, m > 1. (6.2) 
fc=0 k,i=o\ ^'^+^ / 

This equation can be solved perturbatively in g: First we note that for g = 0, the 
model under consideration is nothing but the Gaussian one-matrix model for which 
the function W{z) is known to be 

W^'\z) = j-(^z-^z^-Agy (6.3) 

To the zeroth order in g we hence have 

rTi(o) (2m)! 

m! (m + 1)! 



Using as initial condition ( |6.4D and the normalisation condition Tq = 1, it is obvi- 
ous that equation ( |6.2|) allows one to calculate the individual moments order by or- 
der in g. We remind the reader that just knowing the moment Ti (for n = 1) one 
has the solution of the three-colour problem on a random lattice (cf. equation ( |2.7] )). 
Writing Ti = g^ it appears that to determine Ti'^ one needs to calculate 

{I^^*^}p=i^,.,^j_i;q=i^...^j_p. Below we give Ti to the first six orders in g 

Ti'^ = 2n, 

= lOn + 4n^ 

Tf^ = 70r^ + 60r^2 + 8n^ 

ri^^ = 588n + 764n^ + 240n^ + 16n\ 

Tf) = 5544n + 9520n2 + 4840n^ + 800n^ + 32n^ 

Tf^ = 56628n + 119704^2 + 84216n=^ + 23440n^ + 2400n^ + 64n^ 

Unfortunately, we have not been able to express Ti{g) in a closed form, neither for 
n general, nor for n = 1. However, we can do better than determining the moments 
individually. As we shall see, by an appropriate ansatz we can replace the quadratic 
recursion relation (|6.2| ) by a linear one and solve simultaneously for all Tj. In stead 
of solving the model iteratively in g we shall solve it iteratively in where a is the 
endpoint of the cut of W{z). It is obvious from the eigenvalue integral ( p.2|) as well as 
from the relation (|6.3| ) that small values of g correspond to small values of a^. 

Let us now set t2m = %nd"^ where % = 1 and d = . Then we get from (|6.1|) 
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where 

X = ^ = ^. (6.6) 
9 "^9 



Equation (|6.5| ) must be supplemented by the following boundary condition 

%n 

— ~ const, as m — >• oo. (6.7) 

This relation expresses the fact that W{z) has a square root branch point ai z = a. 
Equation ( |6.5|) can be solved iteratively in d. To leading order in d the solution is again 
given by ( |6.3[) , i.e. we have 



ri'^ = r™, = 1. (6.8) 

Based on an analysis of the structure of the solution found after the first steps of the 
iteration process we introduce the following ansatz 



oo 



%n = '^rn+j Vj d^ , (6.9) 

j=0 

oo oo 

X = Y.X^'^d\ v,=Y.vfd' (6.10) 

j=0 i=0 

with 

= 1, vf^ =0, s ^ 0. (6.11) 

This ansatz fulfills the boundary condition ( |6.7D to any finite order in d and the re- 
quirement that equation ( |6.5| ) should be satisfied and that % = 1 will give us a linear 
recursion relation for the fj's. We note that the function X = X[d) encodes the in- 
formation about the critical behaviour of the model. Critical behaviour occurs when 
the radius of convergence of the power series X = X{d) is reached. We know that we 
will encounter one type of singularity when d approaches ^ but there might be others. 
Inserting the ansatz ( |0|) in (|0| ) we get 



oo 



s=0 [ q=0 q=0 1=0 

-2r^^EE^W^^/+^^'^'(%"^ + ^^)| = 0. (6.12) 

Since this equation must be fulfilled for any value of m and to any order in d the 
quantity in the curly bracket must vanish. This implies that for all s we have 



Y Vs-q \ 5g,0 X - Vq + 2Y Vq+l+1 U d 
q=0 [ 1=0 



l+l 



2nrff^rf'+^,r,+, (^2g + 2/+l^| ^ ^ ^^^^^^ 
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Since vo this equation can only be satisfied for all s if the term in the curly bracket 
vanishes for all q, i.e. 

6g,oX - t;, + 2 £ v,+i+, Ti rf'+i - 2nd £ v, r;+, d'^' ( + 1 ^ ^g^^^^ 

1=0 b,l=0 \ ^(1+ ^ ) 



Now we have succeeded in replacing the quadratic recursion relation ( 5.11 ) by a linear 
one! The relation ( |b.l4 ) must be supplemented by the normalisation condition 



To = ^ TiVid'' = 1. 



(6.15) 



i=0 



From (|6.14| ) and ( |6.15|) we can determine X, {vi}^i and hence W{z) to (in principle) 
any order in d. Inserting the expansion ( |6.1CI| ) in ( |6.14[ ) and ( |6.15| ) we get 

i-1 



(6.16) 



q,0 



,(0 



,(i) 



i-1 

i-1 i-i-1 
Z=0 p=0 



(p) 



A-p-l 'i-p-1 



2g + 2/ + 1 
2g + l 

(0) 



(6.17) 



where on our way we have made use of the initial condition v\ =Q for i 7^ 0. To solve 



recursively these equations we determine in each step first , next X*^*^ and then 
for j > 0. We note that in order to determine X{d) and W{z) to order s in ci we 
need only to know {fp''}i=o,...s;j=o, Since by our iteration process we find indeed 



a perturbative solution of equation (6^) and since this equation is known to have a 
unique perturbative solution our ansatz is justified. Below we give the function X{d) 
to the first seven orders in d 



2n, 

12n, 

lOOn, 

980n - 12n^ 

10584n - 360n^ 

121968n - 7424n^ 

147232872 - 132160^2 + 240n^ 



We remind the reader that the function X = X{d) is related to the endpoint of the 
cut of W{z) by the relation (|6.6| ) and that the three-colour problem corresponds to the 
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case n = 1. Unfortunately, we have not been able to express X{d) in a closed form, 
neither for n general, nor for n = 1. However, we recognise the coefficient of the term 



linear in n as the integral of an elliptic function |11] 

°° \(9kV]'^ 1 /-led 

2rf+12d^+100rf^+980c^'+... = V r/fe = --i-/ K{t^'^)dt-l (6.18) 

(A;!)3(A; + 1)! Svrrf 7o 

We here explicitly see that the model becomes singular when d approaches as argued 
earlier. By means of our results for X{d) we can determine the corresponding value of 
the coupling constant, namely 

- = -J^' 4=-. (6.19) 
gc dc 16 

In particular, using ( 6.18|) , we find to the first order in n 



- = 16\l + n(- -l) +0{n^)\ . (6.20) 



TT 

7 Conclusion 

Exploiting the matrix model formulation of the three-colour problem on a random 
lattice we have developed an algorithm which allows us to solve the problem recursively. 
In addition we have exposed the analyticity structure of the problem and argued why 
the full three-colour problem is so much more difficult to solve than its restricted 
counterpart. 

Let us finish by mentioning an interesting observation. By arguments analogous 
to those given in the introduction for the three-colour model it can be seen that the 
model given by ( p^.3| ) is equivalent to an 0{2n) model on a random lattice where the 
loops are restricted to having even length. We have seen that the critical behaviour 
of the model (|2.3| ) is the same as that of the 0{n) model on a random lattice. Hence, 
restricting the loops to be of even length takes the 0{n) model to the universality class 
of the 0(f) model. 
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